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Engineering complex non-Abelian anyon models with simple physical systems is crucial for topo¬ 
logical quantum computation. Unfortunately, the simplest systems are typically restricted to Majo- 
rana zero modes (Ising anyons). Here we go beyond this barrier, showing that the Z 4 parafermion 
model of non-Abelian anyons can be realized on a qubit lattice. Our system additionally contains 
the Abelian D(Za) anyons as low-energetic excitations. We show that braiding of these parafermions 
with each other and with the D{fLf) anyons allows the entire d = 4 Clifford group to be generated. 
The error correction problem for our model is also studied in detail, guaranteeing fault-tolerance of 
the topological operations. Crucially, since the non-Abelian anyons are engineered through defect 
lines rather than as excitations, non-Abelian error correction is not required. Instead the error 
correction problem is performed on the underlying Abelian model, allowing high noise thresholds to 
be realized. 


I. INTRODUCTION 

Non-Abelian anyons exhibit exotic physics that would 
make them an ideal basis for topological quantum com¬ 
putation m- It has recently become apparent that 
truly scalable quantum computation with non-Abelian 
anyons can only be achieved when invoking active er¬ 
ror correction, despite the protection provided by a finite 
anyon gap HU. The development of practical systems 
in which non-Abelian anyons may be created, manipu¬ 
lated, and detected is therefore highly important. Sys¬ 
tems in which non-Abelian anyons arise typically suffer 
from one of two drawbacks: either they are experimen¬ 
tally extremely challenging to realize (as is the case for 
quantum double [I] or string-net models 0) , or it is not 
clear how they can be made compatible with the active 
error correction required for fault-tolerance (as is the case 
for FQH systems). 

A particularly attractive approach for building a fault- 
tolerant quantum computer is to use a system of physical 
qubits (spin-i particles). A number of technologies al¬ 
low for precise qubit control, such as superconducting 
qubits [8], trapped atomic ions [5], spin qubits uni, or 
cold atoms or polar molecules in optical lattices PH- A 
qubit lattice with two-body nearest neighbour interac¬ 
tions would therefore be an ideal system to realize non- 
Abelian anyons. While the model we propose involves 
three-body interactions, we discuss how these can be ob¬ 
tained from two-body interactions. 

Non-Abelian anyons supported by a qubit system typ¬ 
ically are Majorana zero modes, also known as Ising 
anyons mm- a variety of proposals for experimen¬ 
tal realization of Majorana zero modes in solid state sys¬ 
tems have also been developed m- These anyons can be 
used to perform universal quantum computation when 
enhanced by non-topological operations [17, HHJ- How¬ 
ever, these additional operations are highly resource in¬ 
tensive. Anyon models with a richer set of topological op¬ 
erations would therefore be much more practical for the 
realization of topological quantum computation. Here we 
solve this by introducing a model composed of two-qubit 


Hamiltonian interactions that can realize a more complex 
model of non-Abelian anyons, known as Z 4 parafermions. 
The error correction problem for these is studied in de¬ 
tail. 

Parafermion modes are generalizations of Majoranas 
whose fusion and braiding behavior is more complex and 
computationally more powerful. This has led to a quest 
in recent years for systems that could host them. Nu¬ 
merous proposals for their experimental implementation 
in condensed matter systems such as fractional quantum 
Hall systems, nanowires, or topological insulators have 
recently appeared [TMMj . 

Extrinsic defects in Abelian topological states can be¬ 
have like non-Abelian anyons |35j . The idea of non- 
Abelian anyons at the ends of defect lines, first intro¬ 
duced for FQH states [31j, has been adapted to the D(fLd) 
quantum double models in Refs. isa Eg. These anyons 
are Majorana zero-modes for d = 2 and more powerful 
parafermions for d > 2. Unfortunately, the generalized 
Pauli operators appearing in the D(Zd) quantum dou¬ 
ble models models coincide with the physically relevant 
spin-operators only for d = 2. Otherwise, their struc¬ 
ture makes them highly difficult to realize experimentally. 
The case d = 4, however, allows us to combine the best 
of both worlds. The joint Hilbert space of two qubits 
allows the 4-dimensional generalized Pauli operators to 
be expressed in terms of two-qubit operators. Using this, 
we show how Z 4 parafermions can emerge in a lattice of 
qubits with nearest-neighbor interactions only. This al¬ 
lows the computational power of Z 4 parafermions to be 
harnessed in a qubit system. 

The fact that our system is built on top of a system 
supporting Abelian anyons (the D( Z 4 ) quantum double 
model) proves very useful. The non-Abelian parafermion 
modes can not only be braided with each other, but also 
with Abelian excitations of the quantum double model, 
allowing us to generate the entire Clifford group for d = 4 
by quasi-particle braiding. This extends beyond the lim¬ 
ited set of gates found using the same parafermions in 
previous work |21j . Furthermore, we do not have to per¬ 
form non-Abelian error correction (a still poorly under- 
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stood problem BIEIEM) to guarantee fault-tolerance, 
but can correct the underlying Abelian model. This 
Abelian error correction problem is nevertheless more in¬ 
volved than the well-studied error correction problem for 
the standard D(Zd) models, and we study it in detail. 

The rest of this paper is organized as follows. In Sec.|TT| 
we show how Z\ parafermion operators can be expressed 
in terms of qubit operators. Sec. |III| introduces a qubit 
Hamiltonian whose low-energetic excitations correspond 
to the D(Z 4 ) quantum double model. In Sec. IV we dis¬ 
cuss how Z 4 parafermion modes appear at the ends of de¬ 
fect strings in our model. We demonstrate in Sec. [V] how 
the non-Abelian braiding statistics of these modes can 
be used to perform logical gates. Appendix [D] contains a 
proof that the set of gates which can be performed this 
way generates the entire Clifford group C 4 , which may be 
of independent interest. In Sec. |VI| we study the error 
correction problem of our model in detail and conclude 
in Sec. IVII1 


II. Z 4 PARAFERMION OPERATORS IN TERMS 
OF QUBIT OPERATORS 



FIG. 1. Two qubits are located at each vertex of a Kagome 
lattice. Each pair of qubits hosts two Z 4 parafermions. To 
unlock their potential for non-Abelian braiding, two such 
parafermions need to become unpaired, which is achieved by 
adding a defect line to the lattice. These are strings of strong 
local operators acting on qubit pairs (encircled). They cre¬ 
ate unpaired parafermion modes located at their ends (light 
pentagon-shaped regions consisting of a hexagon and a trian¬ 
gle)- 


III. MODEL 


We consider d-dimensional generalizations of the Pauli 
matrices X and Z. These are unitary operators satisfying 
X d = Z d = 1 and ZX = coXZ, where w = e 2 ™! d with 
integer d > 1. If we define Y = u>^ d+1 ^ 2 X^ Z^, we also 
have Y d = 1, XY = uiYX , and YZ = uZY. Operators 
Xj and Zi act on qudit i and hence [W, Xj] = [ Zi , Zj\ = 
[X,.Z, - ll if, f j. 

These operators are related to those of parafermions. 
Given a total ordering on the qudits {*}, one can obtain 
parafermion operators via a non-local transformation |43| 

72,-1 = (Q X j)Zi, 72, =w (d+1)/2 (n Xj)Zi. (1) 

j<i j<i 

These satisfy the Zd parafermion relations, 

7,'=1, 7*7* = w sgn(fc - j) 7fc7u (2) 

The operators X , Y , and Z can be represented as d- 
dimensional matrices. It is thus natural to seek a rep¬ 
resentation of these operators for the case d = 4 on the 
Hilbert space of two qubits (spins- 5 ). Indeed, given two 
qubits 1 and 2 , one easily verifies that the operators 

X = \( a i + °2 - + ia l a 2) 

Y = i e l37r / 4 (crf + + zcrfcrf + erf erf) 

Z = - ia* 2 ) (3) 

are 4-dimensional generalized Pauli operators, and Z 4 
parafermions can be obtained from these via Eq. 0 . We 
also note that X 2 = erf erf, Y 2 = erf erf, and Z 2 = af erf. 


We consider a two-dimensional Kagome (trihexagonal) 
lattice as in Fig. [T| Each vertex of the lattice hosts one 
4-dimensional qudit (one pair of Z 4 parafermions) or, in 
other words, two qubits. The Hamiltonian of our model 
is given by 

H = Y J H^ + hY J W[+<&)■ ( 4 ) 

A i 

Here, the first term is a sum of equivalent terms for each 
triangle in the Kagome lattice. We label the vertices 
around one triangle a, b , and c, and the two qubits which 
are present at vertex a are called a± and 02 , etc. The 
triangle terms in the Hamiltonian are then given by 

Ha = + a> b >* 2 ) 

- Wal«2 + + <2^cl 

+ <2^cl + <2^ b l<2 + <W b Vc 2 ) ■ (5) 

The second sum ]T7 in Eq. ([4]) runs over all vertices in 
the lattice. The two qubits located at vertex i are called 
il and i2. This second sum thus represents a uniform 
magnetic field in ar-direction. 

Our Hamiltonian involves three-qubit terms of the 
form It is in principle straightforward to gener¬ 

ate these from one-body terms and two-body interactions 
by use of perturbative gadgets [iil - BH] . Consider a “me¬ 
diator qubit” u coupled to qubits a, b, and c. Starting 
from a Hamiltonian 

-Hgadget = —7 2 a u + a ( a a + a b) a u + P a c a u + 7 °"a a b + & a c i 

( 6 ) 
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and consider the perturbative regime A |a|, \ j3\. In 
this regime, it is possible to integrate out qubit u. Tak¬ 
ing up to third-order terms into account, one finds an 
effective Hamiltonian 

H eS = 03 + S)<r* c + (-2^ + 7Kag - 4 °^-* z a a z b a* . 

(7) 
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Choosing 6 = —(3 and 7 = 2^- produces the desired 
three-qubit term without any undesired one- or two-qubit 
terms. 

The generation of three-body interactions in optical 
lattices has been discussed in detail in Refs. 03 Eg. 
These proposals would make the perturbative gadgets 
unnecessary. A “toolbox” for generating spin-lattice 
models such as ours in optical lattices has also been de¬ 
veloped 0H]. Generating non-Abelian anyons other than 
Majorana zero modes by use of perturbative gadgets from 
two-body interactions has previously been discussed in 

Refs, mm- 

The spin-Hamiltonian in Eq. Q can be exactly rewrit¬ 
ten as 

H = -jJ2(Z a Z b Z c + H.c.) + hJ2{Xi + Xj ). ( 8 ) 

A i 

Here again the first sum runs over all triangles in the 
lattice and the corners of a triangle are labeled a, b , and 
c. The second sum runs again over all vertices of the 
lattice. 

We now consider the perturbative limit h <C J and 
regard the second sum in Eq. (| 8 | as a perturbation to 
the first term. Note that all terms in the first sum in 
Eq. ([ 8 ]) commute, so the unperturbed Hamiltonian is 
trivially solved. The lowest-order non-vanishing terms 
appear in sixth-order perturbation theory. We find an 
effective Hamiltonian 

H e ff = - J J2(Z a Z b Z c + H.c.) 

A 

6 Q u6 

- y Qjy £( Xr x l x tXlX v Xi + H.c.), (9) 

where the second sum runs over all hexagons in the 
Kagorne lattice and r, s, t, u, v, w label the six ver¬ 
tices around each hexagon. The effective Hamiltonian in 
Eq. (|9j) is derived in Appendix [A] 

We note that all summands in H e s commute, so the 
system is exactly solvable. The excitations of this system 
are Abelian anyons corresponding to the D( Z 4 ) quantum 
double model. The topological degeneracy of the model 
can be made manifest by studying non-local loop degrees 
of freedom that commute with all stabilizers Z a Z b Z cl 
X r XjX t XlX v X'lj , and their Hermitian conjugates, and 
fulfill themselves Z 4 relations. A possible choice of oper¬ 
ators is illustrated in Fig. [2] 

In passing, we note that the Z 2 version of Eq. ([ 8 ]), in 
which the Z 4 operators X and Z are replaced by Pauli 



FIG. 2. Two sets of logical operators X\ = XX^XX^ ..., 
Z x = ZZZZ.. . (left figure) and X 2 = XX ] XX ] . .., Z 2 = 
ZZ^ ZZ^ ... (right figure) that satisfy the commutation rela¬ 
tions of 4-dimensional generalized Pauli operators. 

operators a x and cr z , leads to an effective Hamiltonian 
analogous to Eq. (|9| and thus provides a very simple 
model with topological order. While this model requires 
three-body operators a z a z a z as opposed to Kitaev’s hon¬ 
eycomb Hamiltonian 1121 which involves two-body inter¬ 
actions only, all of these interactions connect the same 
spin-component, which may provide a significant practi¬ 
cal simplification over the honeycomb model. 

IV. PARAFERMION MODES AND DEFECT 
LINES 

The model is constructed from the cyclic qudit oper¬ 
ators Z and X, which are related to parafermion opera¬ 
tors. It is therefore natural to seek an interpretation of 
the model in terms of parafermionic modes. 

To do this we must first fix the exact form of the 
stabilizers, which define the anyonic charge carried by 
each excitation. Let us use E p (M p ) to denote the sta¬ 
bilizer for a hexagonal (triangular) plaquette, p. For 
hexagonal plaquettes we use the convention that E p = 
XrXlXtX^XyXl,, where r refers to the top-right corner 
and the other corners are labelled in an anti-clockwise 
fashion. For triangular plaquettes we use M p = Z a Z b Z c 
for all triangles of the form A and M p = Z\Z' b Z\ for 
all triangles of the form 77 • The stabilizer operators 
E p and M p are unitary operators with eigenvalues 
k £ {0,1,2, 3}, where here and in the following u> = i 
for d = 4. An eigenvalue u> g of the Ep corresponds to a 
charge anyon of the form e g , while Mp similarly detects 
flux anyons rrih ■ Fusion of charge anyons forms a repre¬ 
sentation of Z 4 , as does that of fluxes. The convention 
for the stabilizer operators chosen before ensures that 
the anyonic charge of both charge and flux type anyons 
is independently conserved (modulo 4). A full clockwise 
monodromy of an e g around an rrih, or vice versa , yields 
a phase u; 9h , see Fig. [ 3 ] for illustration. 

Just as Majorana modes (Ising anyons) in the qubit 
toric code [14, [38], parafermions appear in our system 
at the ends of defect strings. For the interpretation in 
terms of parafermions, it will be useful to introduce a 
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FIG. 3. Phases obtained by braiding the e- and m-excitations 
of the D( Z4) model around each other (top), and by braiding 
the excitations ip and r of the transformed stabilizers around 
each other (bottom). 


new set of composite anyons defined as ip g = e g x m g . 
These also obey Z 4 fusion with each other, and their 
braiding behavior can be inferred from the behavior of 
the constituent charge and flux particles. The particles 
{ipo,ipi,ip2,'4 > 3} form a chiral Abelian anyon model with 
Chern number v = 2 '12!. 

Note that 


e g x m h = ip g x m h - g . 


( 10 ) 


We now perform a local transformation from the set of 
stabilizer generators {E p ,M p }, detecting the charges on 
the left-hand-side of Eq. (10), to a new set {S p , R p } which 


detects the two charges onthe right-hand-side. Let H de¬ 
note the set of hexagonal plaquettes and T denote the set 
of triangular plaquettes. Note that \T\ = 2\H\. Consider 
an injective map p : H —> T, which to each hexagonal 
operator E p assigns one of the six adjacent triangular 
operators M, f ( p y Typically, we choose M p ( p ) to be the 
top-right neighbor of E p , while other choices become nec¬ 
essary next to defect lines. The transformation from the 
old to the new set of stabilizers reads S p = E p for p £ H 
and 


R p = 


M, 


r p E l-H P ) tfpelm(p) 

M p if p ^ Im(yj) 


( 11 ) 


for p £ T. Here, Im(<p) denotes the image of the map ip. 

Since ripen = IIpeT^p = 1, the charges detected 
by the new stabilizers are separately conserved (modulo 
4). Just like the ip g anyons detected by the S p stabilizers, 
the R g charges detected by the R p stabilizers also form 
an anyon model obeying Z 4 fusion. However, while these 
two anyon models have the same fusion rules, they are 
not equivalent, as they exhibit different braiding behav¬ 
ior. A full clockwise monodromy of a ip g around a ip^ 
gives a phase of u } 2gh , a monodromy of a r g around an 
gives a phase of 1 , and a monodromy of a ip g around a 
rh gives a phase of u 9h , see again Fig. [3j Just like the e g 
and m.h charges, the ip g and rh particles correspond to a 
way of decomposing the D(Z 4 ) model into two submod¬ 
els which are closed under fusion, but have non-trivial 
mutual braiding behavior, 


D{ Z 4 ) = {e 0 ,ei,e 2 ,e 3 } x {m 0 ,mi,m 2 ,m 3 } 

= {00,01,02,03} x {ro,ri,r-2,r 3 } , (12) 


where the three particle models other than 
{0o,0i,0 2 ,0 3 } correspond to the simple Z 4 model. 

The stabilizer operators S p detect the presence of ip g 
anyon which are pinned to a pentagon-shaped double 
plaquette, made up of a neighbouring pair of triangu¬ 
lar and hexagonal plaquettes. These anyons can be re¬ 
garded as generalizations of Dirac fermions to the group 
Z 4 (rather than Z 2 ). Just as Dirac modes can be de¬ 
composed into two Majorana modes, so too can the ip 
modes be decomposed into two parafermion modes. Two 
parafermion modes, P a and P&, are therefore associated 
with each double plaquette, P. These are described us¬ 
ing parafermion operators satisfying Eq. Q. The parity 
operator for the ip mode associated with a pair (j. k) is 
defined t f or j <• /., an d so g p = oj 5 ^ 2 'y Pa r yp b . 

For a stabilizer state, the system is within a defi¬ 
nite eigenstate of all Sp. The parafermion modes are 
therefore all paired, with the pairs corresponding to the 
two within each double plaquette. In order to use the 
parafermion modes as non-Abelian anyons, some must be 
allowed to become unpaired. The creation and transport 
of unpaired parafermion modes can be done by adapt¬ 
ing the method of Ref. [35] to the Kagome lattice. The 
method can be interpreted in terms of anyonic state tele¬ 
portation (52) [53], as explained for the Majorana case in 
Ref. [51]. 

The method introduces unpaired parafermion modes 
at the endpoints of defect lines. These are lines on which 
additional single qudit terms are added to the Hamilto¬ 
nian, of one of the two following forms 


+ o-f + of cf + 02) 


Y + H.c. 

lu 5 / 2 XZ^ + H.c. = —^=(of — cr| + of of — ofof ). (13) 




Specific examples are shown in Fig. [4] 

The single qudit terms added along defect lines are 
much stronger than any other interactions, and thus ef¬ 
fectively remove the qudits on which they act from the 
code. This means that the Ep and Mp operators for 
the double plaquettes along these lines no longer com¬ 
mute with the Hamiltonian, and so can no longer be used 
as stabilizer generators. Their pentagon-shaped prod¬ 
uct, Rp , is used instead. The pentagons in Fig. [4] show 
how next to a defect line the mapping ip needs to pick 
the bottom-left triangular-shaped stabilizer of a hexagon¬ 
shaped stabilizer to ensure that their product still com¬ 
mutes with the Hamiltonian. 

This change of the stabilizer generators of the code has 
a drastic effect. Consider an e g anyon moved towards a 
point along a defect line from one direction, and an m g 
moved towards the same point from the other direction. 
Both of these are detected by Rp type stabilizers. When 
they meet on the same double plaquette, they will fuse to 
form a tp g , and so not be detected by the Rp stabilizers 
anymore. In fact, since the Sp stabilizer is removed for 
double plaquettes along a defect line, they will not be 
detected by any stabilizer operator. The ip g occupancy 
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FIG. 4. Strings of alternating single-qudit operators of the 
form w 5/2 Y; + H.c. or oo 5/2 X t Zj + H.c. (encircled) are added 
to the Hamiltonian. These effectively eliminate the qudits on 
which they act from the code, leading to enlarged, pentagon¬ 
shaped stabilizers along the defect string. Parafermion modes 
reside on the pentagons at the ends of the defect strings 
(shaded). A possible choice for two logical operators Xl (top) 
and Zl (bottom) satisfying ZlXl — uiXlZl is illustrated. 


of a defect line corresponds to an increased groundstate 
degeneracy of the system, referred to as a synthetic topo¬ 
logical degeneracy [' 55 ) . 

In the following section, the {'0g ) r /i} decomposition of 
the D{Jj 4 ) model will prove more convenient than the 
{e g , ruh} decomposition. A process in which a defect line 
converts an m g into an e_ g can equivalently be described 
as one in which a r g passes a defect line which emits a 
i>-g- 


V. PARAFERMIONS AS NON-ABELIAN 
ANYONS 


Since unpaired parafermion modes reside at the end¬ 
points of defect strings, it is natural to use them 
to explain the properties of the modified stabilizer. 
Parafermion modes are described by a non-Abelian anyon 
model with particle species {ipo, ipi, V , 2 1 1P3, &}■ Here 
ip o = l corresponds to the anyonic vacuum and a is 
an unpaired parafermion mode. The fusion rules of this 
anyon model are 


a x a = i(>o + ipi + 1P2 + , 

ipg x a = cr, 


(14) 


where ® denotes addition modulo 4. A pair of 
parafermions (or the defect line between them) may 
therefore collectively hold any of the four types of ip 
anyon. 


The anyon model with the fusion rules given in Eq. (14) 


does not allow for a non-trivial solution of the pen¬ 
tagon and hexagon equations. As such, it obeys only 
projective non-Abelian statistics. The computational 
power of braiding parafermions has recently been stud¬ 
ied in Ref. m ■ In our setup, we cannot only braid 
the parafermions with each other, but can also braid the 
Abelian e- and m-particles around them, which provides 
the possibility to perform additional gates. In the follow¬ 
ing, we want to study the gate set that can be generated 
this way. 

As in the Majorana/Ising case, we use four parafermion 
modes (two defect strings) for which the total fusion sec¬ 
tor is vacuum to store one logical qudit. The natural 
logical operators are parity operators for the pairs of 
parafermions. An eigenvalue ui 9 corresponds to a ip g oc¬ 
cupancy for the pair, and so the specific result a x a = ip g 
if they would be fused. We associate the Z basis of the 
logical qudit with the ip occupancy of vertical pairs (con¬ 
nected by defect lines). 

Specific choices of logical operator are illustrated in 
Fig. i The Zl corresponds to a clockwise loop of an e\ 
around a defect line. The braiding of this e\ around the 
ip g held in the pair yields the required phase of lo 9 . The 
Xf corresponds to clockwise loop of an e_i anyon which 
is converted to an m\ through one defect line and back to 
an e_i through the other. Equivalently, we can describe 
it as a clockwise loop of a rr and a transfer of a ipi from 
the right to the left defect line. 

Let us denote a state in which the left defect line holds 
a mode ip g and the right defect line holds a mode iph by 
\tp g ,iph}. Two defect lines create a 4 x 4-fold synthetic 
topological degeneracy. For computational purposes, we 
restrict to the 4-dimensional subspace of states of the 
form \g) L = \ip g , ip~ g ). This is the set of states which can 
locally be created from the anyonic vacuum. The effect of 
the logical operators on these states is Xl \g ) l = |<7 ® 1) L 
and Z L | g) L = u 9 \g) L . 
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FIG. 5. All generators of the single-qudit Clifford group 
can be performed by braiding quasi-particles. The four cir¬ 
cles correspond to the four parafermion modes which are 
used to store one logical qudit. The left part of the fig¬ 
ure illustrates how to perform the logical operators X and 
Z by braiding the Abelian excitations of the D(Z 4 ) model 
around the parafermions. The right part demonstrates a log¬ 
ical Hadamard gate H, which is performed by braiding the 
parafermion modes themselves. 


In addition to the logical operators Xl and Zl 1 which 
can be performed in our model by braiding the Abelian 
D( Z4) anyons around the parafermion modes (ends of 
defect strings), we can perform further topologically pro¬ 
tected single-qudit and two-qudit gates by braiding the 
parafermion modes themselves. Defect lines used for 
braiding are shown in Appendix [Bj Crucially, braiding 
parafermions allows one to perform an entangling gate 
by topological means, which is in contrast to Majorana 
fermions m- What is more, exploiting the fact that our 
non-Abelian system is built on top of an Abelian D(Z 4 ) 
system allows us to generate the entire 4-level Clifford 
group by braiding quasi-particles, as we discuss in the 
following. 

For the rest of this section, X and Z refer to the logi¬ 
cal operators called Xl and Zl before, respectively. The 
first column in Fig. [ 5 ] illustrates how braiding of D( Z4) 
charges and fluxes can be used to perform logical X and 
Z gates. Whether an e\ or an m 1 anyon is used to per¬ 
form the logical Z gate is irrelevant. 

Consider two ?parafermion modes storing a ip g par¬ 
ticle. A full clockwise monodromy of one parafermion 

around the other can be understood as a monodromy 

2 

of the constituent e g around the m g , yielding an u> 9 
phase. We can thus expect a single exchange of the two 
parafermion modes storing a ip g to yield a square root 
of this phase, such as w® 2 / 2 . This is demonstrated di¬ 
rectly by studying the necessary microscopic operations 
in App. [Cj 

For a logical qudit stored in four parafermion modes, 
let S denote a clockwise exchange of a vertical pair 
of parafermion modes, and T an exchange of a hori¬ 
zontal pair, see Fig. [6] As discussed, we have S = 
XOgW 9 '/ 2 |<?) (g |, while T is diagonal in the logical X ba¬ 


sis. In the logical Z basis, T reads (for d = 4) 

T= l -e-™/*Y j e i ^ 9 - h)2 \g){h\ . (15) 

gh 



FIG. 6 . Generators S and T of all gates that can be performed 
on a qudit stored in the fusion space of four parafermions by 
braiding them. 

Again, in contrast to Majorana fermions, parafermions 
support an entangling gate between two logical qudits 
by braiding operations pH]. The controlled phase-gate 
A is defined by its action on a logical two-qudit basis 
state, A | g,h) = uj 9h \g,h). In our parafermion scheme, 
an entangling gate can be performed by braiding of a 
pair of parafermions from one qudit with a pair from 
the other. Let us consider, for example, the braiding 
of the left vertical pair for both qudits. For an initial 
logical product state | g,h), the process corresponds to 
braiding a ip g clockwise around a iph , which yields a phase 
of u} 2ah . The resulting operation is therefore the squared 
controlled phase-gate A 2 . For d = 2, corresponding to 
the Ising/Majorana case, A 2 = 1, and so this operation is 
trivial. For d > 2, however, it is a non-trivial entangling 
gate, akin to the one proposed in Ref. m- 

Clearly a more powerful entangling gate would be A 
itself. This can be achieved for Z^ parafermions for odd 
d by taking the (d+ l)/2-th power of A 2 . However these 
do not admit the simple decomposition into qubits that 
we have used in defining the model. Fortunately, we can 
make use of the underlying charge and flux anyons to 
realize A despite the even qudit dimension. 

The defect line may be interpreted as a hole for ip type 
anyons: an area in which they may be placed such that 
their state becomes delocalized along the line and they 
are no longer detected by the stabilizers [55H55] . Similar 
holes can also be engineered for the constituent charge 
and flux anyons. A defect line is therefore a special case 
of the combination of a charge and flux hole, in which 
only ip g = e g x m g type anyons may reside rather than 
general e g x mh anyons. Nevertheless, we can consider 
a process in which a defect line is transformed into a 
separate charge and flux hole. Details on these holes 
and the transformations between them can be found in 
Appendix [E] 

When only the charge hole of one qubit is braided 
around the defect line of another, the process for an ini¬ 
tial state | g, h) corresponds to braiding an e g around a 
iph , which would yield the phase uj 9h . The charge and 
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FIG. 7. Performance of a controlled phase-gate. The a) part 
of the figure shows a logical product state | g, h) stored in the 
fusion space of eight parafermions. The defect line storing 
a mode ip g can be split into two defect lines storing DifLf) 
charges e g and m g , respectively. Braiding both endpoints of 
one of these lines clockwise around the defect line storing the 
iph mode, as shown in the b) part, produces a phase ui 9h , as 
required. 


flux holes can then be recombined into a defect line. The 
net effect of the entire process is to apply the controlled 
phase gate A. Such a process is illustrated in Fig. [ 7 ] 

One process which could split the defect line in this 
way is simply to intersect it with two others. One would 
be a line along which charge anyons are hopped by high- 
strength terms. The other would similarly hop flux 
anyons. The stabilizers that detect charges and fluxes, 
respectively, along these lines would then be suppressed. 
By adiabatically removing the defect line which delocal¬ 
izes if modes, its if g anyon occupation would be trans¬ 
ferred to these two lines. The recombination of the defect 
line would be done by the reverse process. 

For a tensor product of d-level systems, the Clifford 
group Cd consists of gates that map tensor products of 
d-level Pauli operators to other such tensor products un¬ 
der conjugation. In Appendix [D] we prove the following 
theorem. 

Theorem. The single-qudit gates S, T, and Z, and 
nearest-neighbor controlled phase-gates A generate the 
entire Clifford group C 4 . 

As an example, H = STS = TST satisfies HXH' { = 
Z and HZH* = X\ so it can be identified with the 
logical Hadamard gate, up to a phase. Indeed, using the 
standard definition 

# =4E^I<7>W , (16) 

v gh 


An X corresponds to a transfer of a ifi from the right 
to the left defect line, accompanied by a clockwise loop 
of a r 1 around a horizontal pair. A Z corresponds to a 
clockwise loop of a r\ around a vertical pair. A w/2 rota¬ 
tion as performed by H thus maps these two operations 
onto each other, up to the fact that we do not perform 
a vertical ifi transfer, as the if occupancy of the vertical 
pair is delocalized along the defect line. 


VI. ERROR CORRECTION 

For any system with a finite energy gap at finite tem¬ 
perature, excitations will appear with a finite density. 
This corresponds to finite length scale on which quantum 
computation can be performed before errors are almost 
certain to appear. This length scale can be increased by 
increasing the gap or lowering the temperature. How¬ 
ever, neither of these methods is truly scalable. Error 
correction is therefore required if scalable quantum com¬ 
putation is to be performed. 

For non-Abelian systems, the first studies of the cor¬ 
responding error correction problem have recently ap¬ 
peared H 0 IMS 2 - Error correction for non-Abelian 
anyons is still poorly understood and its feasibility has 
not been demonstrated for the (realistic) time-continuous 
case. It comes thus very welcome that while our sys¬ 
tem provides the computational power of non-Abelian 
parafermions, its physical excitations still are Abelian 
D( Z 4 ) anyons, and the error correction problem for 
D(Z n ) quantum double models (including the time- 
continuous case) is well-studied [39 | l40 | [59ll6l] . However, 
when correcting these D{ Z 4 ) anyons, we face a number of 
difficulties not considered in previous studies [5H1 155H5I] : 

(i) Our stabilizer operators are products of Z 4 qudit 
operators X, X { , Z, Z\ while an error model is 
realistically expressed in terms of single-qubit op¬ 
erators a x , a y , a z . These do not map eigenstates 
of the stabilizer operators to other eigenstates and 
one single-qubit operator can produce a product of 
up to three qudit operators (see below). 

(ii) We consider quantum information stored in a syn¬ 
thetic topological degeneracy, which involves a de¬ 
fect line allowing anyons to change from one sublat¬ 
tice to the other (stars to hexagons and vice versa). 
We thus cannot decode each sublattice separately, 
as usually done for the toric code and other D(Zd) 
quantum double models, but have to correct both 
of them simultaneously while taking the possibility 
of transferring anyons from one to the other into 
account. 


one verifies that = H. 

One possible implementation of H (up to a phase) is 
a cyclic permutation of the four parafermion modes, as 
in Fig. [5} This can be pictorially understood as follows. 


(iii) Besides simplistic i.i.d. error models (such as de¬ 
polarizing noise), we are particularly interested in 
Hamiltonian protection of a quantum state subject 
to thermal errors. 












^.X X 

°1 , °2 

x , x\ xz 2 , x ] z 2 

y y 

xz , x f z, xz\ x f z f 

^2 2 

CT 1 , a 2 

z, z f 


TABLE I. Conversion from single-qubit Pauli operators to 4- 
dimensional generalized Pauli operators. When a syndrome 
measurement is performed, a Pauli operator is converted to 
each of the generalized Pauli operators in the right-hand col¬ 
umn with equal probability. 

(iv) We do not consider a square lattice, but a trihexag- 
onal one, which makes moving anyons and defining 
their distance more involved. 


A. Error model 

Since our 4-level qudits are composed of two qubits, it 
is natural to consider an error model in terms of single¬ 
qubit operations cr^, a y , and a z . For a qudit hosted in 
two qubits 1 and 2, single-qubit Pauli operators can be 
expressed in terms of Z 4 operators by inverting Eq. [3] 
We find 

of = \x{l - Z 2 ) + H.c., 

C 2 = \x{l + Z 2 ) + H.c., 
a\ = i e i57r/4 y(l + Z 2 ) + H.c., 
a v = I e i37r/4 V(l-Z 2 )-f-H.c., 
al=e~ in/i Z + H.c., 

@2 = e^ /A Z + H.c.. (17) 

If we start from an eigenstate of all stabilizer operators, 
applying single-qubit Pauli operators will generate a su¬ 
perposition of states corresponding to different syndrome 
outcomes. By measuring all stabilizer operators, we can 
project again into a subspace with definite syndrome val¬ 
ues. Each single-qubit Pauli operator thereby translates 
into a product of up to three qudit operators. Table |T] 
summarizes (up to irrelevant phases) into which qudit 
operators a certain single-qubit Pauli operator will trans¬ 
late with equal probability. 

As a first simple error model, which does not involve 
a notion of Hamiltonian protection, we consider depolar¬ 
izing noise. That is, for each qubit of the code we apply 
a Pauli operator with some probability p (the depolar¬ 
ization rate), where each of the three Pauli operators is 
chosen with equal probability. 

More interesting from a physical perspective is a ther¬ 
mal error model. We consider a quantum state stored in 
the degenerate groundstates of the Hamiltonian given in 
Eq. and assume that the system is weakly coupled 
to a heat bath at some temperature T. Following e.g. 
Ref. (S] , we assume that evolving the system according to 


the Metropolis algorithm provides a reasonable approxi¬ 
mation of the thermalization process, since the evolution 
obtained by means of the Metropolis algorithm is local, 
Markovian, and has the thermal state as its unique fixed 
point. 

During our simulation, we proceed as follows. We first 
pick one of the spins-^ of the system at random, then 
pick one of the three single-qubit opertors acting on that 
qubit at random, and convert that to a 4-dimensional 
generalized Pauli operator according to Table |T] We then 
calculate the energy cost A to t of applying that general¬ 
ized Pauli operator (or products thereof). This energy 
cost is of the form 

A to t = toA a + nAo , (18) 

where A^ and Aq are the energy costs of creating 
a single triangle/hexagon-type anyon with charge 1 or 
3 in Eq. ([ 9 ]), respectively. (That is, A a = 2 J and 
Aq = 2 ^(yjjs-) Creating an anyon with charge 2 will 
have an energy cost 2 Aa or 2Ao- The coefficients m and 
n are elements of {0,±2,±4}, depending on the change 
in anyonic charge. The proposed error is then accepted 
with probability min{ 1, e _Atot / fcBT }. The noise model we 
apply is thus the standard classical Metropolis algorithm 
that maps eigenstates of Eq. © to other eigenstates. At 
any given time during our simulation, the system is “clas¬ 
sical” in the sense that it does not involve superpositions 
of different anyon configurations. 

If the proposed error is accepted, we copy the current 
state of the system and try to correct it. If correction 
is successful, we continue our simulation with the un¬ 
corrected version of the system. If correction fails (for 
at least one logical operator), we interpret this as the 
quantum information having survived for a time which is 
given by the number of Metropolis steps divided by the 
number of spins in the system. 

The thermal error model has three relevant energy 
scales I’sT, Aq, and Aa- Since Aq appears in higher- 
order perturbation theory than Aa, we expect A a > 
Ao- Furthermore, effective protection requires fcsT < 
Aq , A a ■ We introduce a parameter A which quan¬ 
tifies the separation of these three energy scales, i.e., 
Aq = A ksT and Aa = A 2 ksT. Very high values of 
A are uninteresting, since they exponentially suppress er¬ 
rors from occurring. 


B. Without defects 

If there are no defect lines present, the anyonic charge 
of both types of anyons is conserved (modulo 4), and 
they can be corrected separately. Various techniques 
have been developed for correcting general D(Z n ) quan¬ 
tum double models 391 f59UBTj . However, correcting the 
D(Z 4 ) case is particularly easy, since we can exploit the 
relation Z 4 /Z 2 ~ Z 2 . Specifically, we can first fuse all 
oddly-charged anyons in pairs. In a second round, the 
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FIG. 8. Error rates pl of the logical operators X] and Z\ 
illustrated in Fig. [2] as a function of the qubit depolarization 
rate p for code sizes L = 20, 28,36,44, 52. Each data point 
represents 10 4 logical errors, such that error bars are negligi¬ 
ble. We recognize a threshold error rate p c ft: 24% for Xi and 
p c « 10% for Z\. 


remaining anyons, which are all of charge 2 , are fused in 
pairs. In order to find these pairings, we use the library 
Blossom V [62] , which is the latest implementation of 
the efficient minimum-weight perfect matching algorithm 
due to Edmonds [53] ■ The weight between two equal-type 
anyons is thereby defined as the minimal number of gen¬ 
eralized Pauli operators that need to be applied to create 
a pair of anyons at the two given locations. 

Fig. m shows our results for the depolarizing noise 
model, i.e., the logical error rates of the the logical op¬ 
erators X\ and Z\ illustrated in Fig. [ 2 ] as a function of 
the depolarization rate p. One clearly recognizes thresh¬ 
old error rates p c ft 24% and p c ft 10%, respectively. 
The equivalent figures for the logical operators Z 2 and 
X '2 look very similar and yield equivalent threshold error 
rates p c . 

These thresholds are best compared with those for an 
equivalent code based on Z 2 anyons, and so with only 
a single qubit on each vertex. For independent bit and 
phase flips, the thresholds for X\ and Z,\ are p c ft 16.4% 
andp c ft 6.7%, respectively [64ll65|. When the hexagonal 
and triangular plaquettes are decoded separately, these 
correspond to thresholds of p c ft 24.6% and p c ft 10.5% 
for depolarizing noise. The similarity of these Z 2 values 
with those of Z 4 is remarkable. This qudit code is there¬ 
fore just as adept at suppressing qubit noise as its qubit 
counterpart. 



FIG. 9. Average lifetimes r of the logical qudit with logical 
operators A'i and Z\ as a function of code size L for A = 
1,2, 3, 4, 5. Each data point represents 10 4 experiments. The 
lifetime is defined as the number of Metropolis steps until the 
first logical operator detects an error, divided by the number 
of spins in the code. 

It is well-known that the finite-temperature lifetime of 
a two-dimensional quantum memory with local interac¬ 
tions only is upper-bounded by a constant independent 
of the system size, see e.g. Ref. [ 66 ]. Fig.[9]shows the life¬ 
time of a logical qudit with logical operators A'i and Z\ 
subject to the thermal error model. We notice lifetimes 
that decrease to an asymptotic value for large L and con¬ 
siderable finite-size tails. These tails correspond to the 
regime in which the breakdown of error correction is not 
due to the density of anyons becoming so high that pair¬ 
ing them becomes ambiguous, but where the breakdown 
is caused by one of the first pairs wandering along a topo¬ 
logically non-trival path around the torus. The smaller 
the system, the longer it takes to produce an anyon pair, 
leading to the observed tails for small enough L and T 
(large enough A). 

C. With defects 

When defect lines as in Fig.[4]are present, the error cor¬ 
rection problem becomes more involved. It is no longer 
possible to correct the two anyon types (hexagons and 
triangles in our case) separately, as is usually done for 
the D{ Z„) models (39 ) (591161] . Instead, error correction 
needs to take the possibility of converting between dif¬ 
ferent anyon types into account. We thus pair all oddly- 
charge anyons of both types in a first round and all re¬ 
maining charge 2 anyons of both types in a second round. 
Pairings can involve anyons which are of equal or of dif¬ 
ferent type. The weight for connecting two anyons is 
defined as the minimal number of generalized Pauli op¬ 
erators needed to create a pair of anyons at their respec¬ 
tive positions from the vacuum. For equal-type anyons, 
this will be an error string that crosses an even number 
of defect lines, while for different-type anyons this will 
be an error string that crosses an odd number of defect 
lines. This can mean, for instance, that connecting two 
equal-type anyons can have a large weight despite them 
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FIG. 11. Average lifetimes r of the logical qudit stored in 
the defect logical operators X and Z illustrated in Fig. [ 4 ] as 
a function of code size L for A = 1, 2,3,4, 5. Each data point 
represents 10 4 experiments. The lifetime is defined as the 
number of Metropolis steps divided by the number of spins in 
the code. 
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VII. CONCLUSIONS 


FIG. 10. Error rates pl of the logical operators Xl and Zl 
illustrated in Fig. [4] as a function of the qubit depolarization 
rate p for code sizes L = 20, 28, 36, 44, 52. Each data point 
represents 10 4 logical errors, such that error bars are negli¬ 
gible. We recognize a threshold error rate p c ss 24% for Xl 
and p c « 10% for Zl- 


being geometrically nearby, if there is a defect line be¬ 
tween them. 

For a code of linear size L in both dimensions, with 
periodic boundary conditions and L even, we choose de¬ 
fect lines involving L/2 + 1 qudits, as shown in Fig. [4] for 
L = 20. The logical operators Xl and Zl then have a 
distance L + 2 and L/2 + 4, respectively. 

For the depolarizing error model, we find the threshold 
error rates p c for both of the logical operators Xl and 
Zl given in Fig. [ 4 ] The results are given in Fig. [To} For 
the defect operator Xl, we find a threshold error rate 
p c ~ 24%, as for the operators X\ and X 2 in the defect- 
free case (Figs. [ 2 ] and [8]), while for the defect operator 
Zl we find a threshold error rate p c « 10%, as for the 
operators Z\ and Z 2 in the defect-free case. 

The fact that these values coincide with the defect free 
case is not unexpected. The introduction of the defects 
essentially corresponds to a change in the boundary con¬ 
ditions. However, the vast majority of errors have large 
support within the bulk. The value of the threshold is 
therefore dominated by bulk effects rather than bound¬ 
ary effects. 

Fig- 0 shows the average lifetime of the qudit stored 
in the synthetic topological degeneracy in Fig. [4] for the 
thermal error model. We note that for a given parameter 
A, the asymptotic lifetimes (L —» 00 ) are close to those 
in the defect-free case given in Fig. [9] 


We have proposed a system which, on the physical 
level, involves only nearest-neighbor two-qubit interac¬ 
tions, allows one to perforin all Clifford gates through 
quasi-particle braiding, and has a well-understood error 
correction problem. 

We have greatly benefitted from the fact that our non- 
Abelian system is built on top of a system whose ex¬ 
citations correspond to an Abelian anyon model. This 
allows us to perform the logical operators X, Z, and A 
through quasi-particle braiding. It also makes our error 
correction problem manageable, despite some subtleties 
such as the fact that single-spin Pauli operators gener¬ 
ate superpositions between different syndrome outcomes 
and the ability to convert between different anyon species 
during error correction. 

Universal quantum computation requires the ability to 
perform non-Clifford gates, such as “small-angle” uni- 
taries. While it is not difficult to perforin a non-Clifford 
operation by non-topological means in our system, this 
abandons fault-tolerance. The technique of magic state- 
distillation j67| is typically used to restore fault-tolerance. 
While research on magic state distillation has so far 
focused on prime qudit dimensions d |68l 159] . qudit 
codes with the right transversality properties to perform 
magic state distillation in non-prime dimensions, includ¬ 
ing d = 4, also exist EH- Unfortunately, for non-prime 
d it is not known whether Clifford gates plus an arbi¬ 
trary non-Clifford gate are sufficient to achieve univer¬ 
sality m- It is our hope that our work fuels interest in 
the d = 4 case, being a power of 2 and thus allowing to 
employ qubits, as demonstrated in our work, while being 
the smallest power of 2 that allows one to go beyond the 
Ising/Majorana case. 

Alternatively, one could imagine energetically penaliz¬ 
ing one of the degrees of freedom of a two-qubit Hilbert 
space to obtain a synthetic qutrit (d = 3). Magic state 
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distillation for qutrits is well-studied m , potentially al¬ 
lowing to perform fault-tolerant universal quantum com¬ 
putation with Z 3 parafermions in a qubit system m- 
The authors thank M. Barkeshli for elaborations on 


the development of the idea of generating non-Abelian 
defects in topological systems. This work was supported 
by the SNF, NCCR QSIT, and IARPA. 
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Appendix A: Sixth-order degenerate perturbation 
theory 

For our perturbation theory, we employ a Schrieffer- 
Wolff transformation [74], as formalized in Ref. [75]. 

Consider an unperturbed Hamiltonian H 0 whose spec¬ 
trum can be separated into a low- and a high-energy sub¬ 
space, which are energetically separated by a gap. Given 
a perturbation V , we want to find an effective Hamilto¬ 
nian H e g describing the “effective” physics on the low- 
energy subspace. The effective Hamiltonian can be de¬ 
veloped in a perturbative series 

H eB = H™+H& ) +H™ + ... (Al) 

in powers of some small expansion parameter. 


Let P denote the projector onto the low-energy sub¬ 
space and Q = 1 — P the projector onto the high- 
energy subspace. We define Vd = PVP + QVQ and 
V>d = V — Vd = PVQ + QVP. For some operator A, 
we define the superoperator A via A(O) = [A, 0\. Let 
Ho = J2iEi\i)(i\ be the spectral decomposition of H 0 
and define the superoperator C via 

C(0) = ]T - H.c.. (A2) 

i,j 1 3 

We employ the convention that unless indicated other¬ 
wise by use of brackets, a superoperator C acts on all 
operators to its right. 

For the sixth-order effective Hamiltonian, one derives 
from Ref. [75] the expression 


H. 


( 6 ) 


eft = ^PS 5 {V od )P - —P( s i s 3 + S&S 1 + S 3 Sf + S^S 1 + S 2 S 1( S 2 + SiSfiiVoJP + —PSj \V od )P , 


(A3) 


where 


51 = C(V, 3d) 

5 2 = -£V d (Si) 

5 3 = -CV d (S 2 ) + ^CSl(V od ) 

5 4 = -CV d (S 3 ) + ^£(SiS 2 + S 2 S 1 )(V r od ) 

= -CV d (S A ) + + SJa + SaSO (V od ) 

- ^£St(Vod) • (A4) 


In our case, the low-energy subspace onto which P 
projects is given by the space in which all triangle op¬ 
erators in Eq. ([8]) have minimal energy, i.e., Z a Z^Z c = 1 
for all triangles (a, b, c). This subspace is fully degener¬ 
ate. The lowest-energetic excitations change the eigen¬ 
value of a stabilizer Z a Zt,Z c from 1 to ±i. Since the 
eigenvalue of — (Z a Z b Z c + H.c.) is thereby changed from 
—2 to 0, this has an energy cost A = 2 J. Note, however, 
that stabilizer eigenvalues can only be changed in pairs, 
such that the gap between the low-energetic (ground- 
state) subspace and the space of excited states is in fact 
given by 2A. 


A crucial property of our Hamiltonian is that there is 
no lower-than-sixth-order perturbation that acts within 
the groundstate space. Therefore, we are only interested 
in terms of the form PV od (V d ) 4 V od P, which allows to 
greatly simplify the effective Hamiltonian. Namely, only 
the first summand in all expressions in Eqs. (A3) and 


(|A4|) is relevant in our case. We find 

(CV d ) 4 (CV od ),V od 


1 


/A 6 ) — -p 

«eff - 


P. 


(A5) 


Using now that in our case VdP = 0, this can be further 
simplified to 


= X -PL (£ (£ (£ (£ (U od ) V d ) V d ) V d ) V d ) V od P 

- ^PV od {CV d f{CV od )P 
=-PV od {CV d )\LV od )P. (A6) 


There are 6! = 720 possibilities for applying the six fac¬ 
tors X r XlX t XlX v Xl, around one hexagon which leads 
the system back to the groundstate. Table [H| lists all pos¬ 
sible routes the excitation energy above the groundstate 
can take, together with their numbers of possibilities. 

In conclusion, we find the sixth-order effective Hamil¬ 
tonian 

Peff = -q^{X r X\X t XlX v Xl + H.c.), (A7) 


where the dimensionless prefactor 


96 

48 

48 

48 

-1- 

H-1- 

32 

64 

64 

64 

192 

24 

72 

4---f- 

-1- 

256 

63 

128 

384 


8 


96 

128 


+ 


96 

128 


(AS) 


is given by the multiplicities in Table |TT] divided by the 
product of all excitation energies (in multiples of A) 
along the virtual process. 
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72 


TABLE II. Possible routes the excitation energy above the 
groundstate can take (left column), together with their re- 
spectiv multiplicities (right column). Note that the number 
of multiplicities adds up to 6! = 720. 



FIG. 12. A selection of double plaquettes used in a 
braiding operation. Each double plaquette corresponds to 
two parafermion modes. For a double plaquette P, the 
parafermion to the right is labelled Pi , and that to the left is 

Pi. 


Appendix B: Moving unpaired parafermion modes 


To consider the creation and braiding of unpaired 
parafermionic modes, we must first decide on the dou¬ 
ble plaquettes with which we will work. Let us consider 
those of Fig. [12] To visualize the two parafermion modes 
within each double plaquette we use light blue circles. 
The one to the right of a double plaquette P is labelled 
Pi, and that to the left is P 2 . 

Parity operators for ip modes are defined on pairs of 
parafermion modes. We are primarily concerned with 
two types of pairing: those of the two modes within the 
same double plaquette, and those of two modes from 
neighbouring double plaquettes. Relevant examples of 
the latter type are shown in Fig. [T2] by red, orange and 
green lines connecting the corresponding modes. 

For the two modes within each double plaquette, the 


parity operator w 5 / 2 7 p 1 7 p 2 corresponds to the stabilizer 
Sp. The orange and red lines connecting modes P 2 to 
(P + l)i denote the parity operators w 5 ^ 2 7p 2 7( P _ t _ 1 ) 1 ■ For 
orange lines, these correspond to the operator Y on the 
vertex through which the line passes. For red lines they 
correspond to the operator X' Z. 

Consider a state initially within the stabilizer space. 
The parity operators for the pairs of parafermion modes 
within each double plaquette are therefore part of the 
stabilizer. The state therefore corresponds to this definite 
pairing of the modes. 

Let us now consider the removal of the operator Sa 
from the set of stabilizer generators (while Ra remains). 
The corresponding parafermion modes are now, in some 
sense, unpaired. This contributes a factor of four to the 
ground space degeneracy m ■ However, due to the fact 
that the ‘unpaired’ parafermions are not well separated, 
it is not difficult for local perturbations to lift the de¬ 
generacy of this space. To become truly unpaired, and 
benefit from topological protection, they must be sepa¬ 
rated. 

To do this, we can add a term K(Y + W) to the 
Hamiltonian, which corresponds to the parity operator 
w 5 / 2 7a 2 7bi- This acts on the vertex through which the 
orange line connecting these modes passes. 

For K J, this new term will overwhelm the Sp term. 
The pairing of Bi and B 2 will then be broken, and Bi 
will become paired with A 2 instead. The unpaired mode 
originally at A 2 is therefore effectively moved to B 2 . If 
the new term is introduced adiabatically, the degener¬ 
ate subspace associated with the unpaired parafermion 
modes will remain in the same state during this process. 

Similar processes can be used to move the unpaired 
modes further. The Hamiltonian term K{X' Z + XAZ) 


corresponding to can then be used to move 

the parafermion at B 2 to C 2 , for example. Unpaired 
parafermion modes can therefore be separated by arbi¬ 
trary distances, at the endpoints of lines on which sin¬ 
gle qudit terms are added to the Hamiltonian. In terms 
of qubits, these correspond to two-body interactions be¬ 
tween qubits in the same site. 

In order to unlock the potential of parafermions for 
quantum computation, it must be possible to braid the 
parafermion modes. Let us consider a specific example 
of this, using the system of Fig. 12 Consider an initial 
state within the stabilizer space of all Sp except A and L. 
At these two points, we have the unpaired parafermion 
modes A±, A 2 , L\ and L 2 . Let us now consider operations 
such as those described above to move A\ and L 2 away, 
beyond the bottom of the figure. All four parafermion 
modes are then well separated, and so the ground state 
degeneracy is topological protected. 

We will now consider the exchange of A 2 with L\. We 
do this by first moving A 2 to K 2 , then L\ to A 2 , and 
finally K 2 to L\. The two modes have then swapped 
places. An exchange of opposite chirality would corre¬ 
spond to first moving Li to Ki, and so on. Note that all 
modes are kept well separated during the exchange, and 
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FIG. 13. Double plaquettes used in the worked example of a 
single exchange. 


so topological protection is always maintained. 

During the exchange, the movement of the modes is 
mostly achieved using Hamiltonian terms that corre¬ 
spond to the red and orange pairings in Fig. [12] These are 
all single qudit terms. At the junction, however, terms 
corresponding to the green pairings are used. We must 
therefore consider these in detail. 

For the pairing show by the light green line, the parity 
operator is w 5 ,/ 2 7 p 2 7 q i . This has the effect of creating 
an ip g , ip_ g pair on the double plaquettes O and P. This 
requires the two qudit operator X\X 4 Z\. For the dark 
green pairing, the w 5 / 2 7 p 1 7 g . 9 parity operator similarly 

requires the three qudit operator w 5 / 2 Xl x\Z 2 Zz. These 
terms correspond to four- and six-body quasi-local inter¬ 
actions on the corresponding qubits, respectively. They 
can be realized by standard methods of perturbative gad¬ 
gets. However, note that they need only be implemented 
while an exchange is in progress. 


Appendix C: Exchange of two parafermion modes 

To determine the effects of a single clockwise exchange 
we consider the smallest possible implementation. This 
involves the double plaquettes labelled F, G and P in 
Fig. [T2] which are shown in more detail in Fig. [13] We 
consider a state in which the modes at P\ and F 2 are un¬ 
paired, and those at Gi and G 2 are paired by the Hamil¬ 
tonian term So- Using this, we determine the effects of 
exchanging the unpaired modes. The method used to 
exchange the two modes is similar to previous methods 
proposed in order to perform anyon braiding [521153] . 

The results of the exchange are most easily understood 
in terms of the ip mode formed by this pair. The parity 
operator, T for this mode is an operator that creates a 
ipi, ip -1 pair and places them in double plaquettes P and 
F, respectively. Also it must have eigenvalues of the form 
uj 9 , and so T 4 = 1. These conditions are satisfied by 

r = uj il+2a)/2 z 1 xlz 2 z 3 . (ci) 

We similarly require operations that can move 
parafermions between the relevant plaquettes. These cor¬ 
respond to the green line between G and P and the red 


line between F and G. These are 

n = ^ 1+2b)/2 XiX\z 2 Z 3 , <F = w ( 1 + 2 c)/ 2 X 1 t Zi, (C 2 ) 

respectively. In these relations a, b and c are all elements 
of Z 4 . 

We have freedom in choosing the values o, 6 , c £ Z 4 
for these relations. The corresponding freedom also ex¬ 
ists for all operators used to move parafermion modes, 
as well as the logical operators. The values used do not 
simply correspond to differences in a global phase. In¬ 
stead they determine which eigenspace of these operators 
has eigenvalue co° = 1 , and so which one corresponds to 
the vacuum occupancy ipo of the ip mode. These phases 
therefore cannot be chosen entirely arbitrarily, since the 
overall conservation constraint of ip modes must be main¬ 
tained. However, since here we do not explicitly consider 
the operations that placed unpaired parafermion modes 
at Pi and F 2 , we can assume that their phases are cho¬ 
sen in a way that maintains this conservation. We will 
therefore consider a free choice of a, 6 , and c. 

The first step in exchanging the parafermions is to 
move the one at Pi to G 2 . This is done by adiabati- 
cally changing the Hamitonian to one in which the term 
n + nl is present and stronger than So- This causes the 
modes at Gi and Pi to pair, moving the mode once at Pi 
to G 2 . The mode at F 2 is then moved to that at Pi by 
adiabatically changing the Hamitonian to one in which 
the + term is present and stronger than Sq, and the 
n + nt term is removed, pairing F 2 with G The mode 
at G 2 is then moved to F 2 by adiabatically removing the 
n + nf term and so allowing Sq to become dominant 
and Gi and G 2 to pair. This process then results in the 
clockwise exchange of the modes. 

The first step of this transformation takes a state that 
is initially in the w° eigenspace of Sq and projects it to 
one in the oj° eigenspace of n. The next step projects the 
state into the w° eigenspace of <f>. The final step projects 
back into the w° eigenspace of Sq- The end effect is 
then PgP^PuPg- Here Pq is the projector onto the w° 
eigenspace of Sg , etc. The rightmost Pq simply reflects 
the fact that the initial state lies within this eigenspace. 

The parity operator F for the pair of unpaired modes 
commutes with Sq, and so can be mutually diagonalized 
with the above operator. We can therefore interpret its 
effects in terms of the phase factor assigned to each of 
the possible ip g eigenspaces of the ip mode of the pair. 

When doing this, different values of a, 6 , and c will 
result in different operations. This may seem to con¬ 
tradict the standard notion of a topological protected 
operation. However, these differences can be most eas¬ 
ily understood by considering movement of parafermion 
modes implemented by measurement rather than adi¬ 
abatic Hamiltonian manipulation. This method forces 
pairing of parafermion modes by measuring the occu¬ 
pancy of their corresponding ip mode, and so forcing it to 
have a definite value. Ideally, this measurement will give 
the vacuum result ipo . The effect is then the same as the 
adiabatic manipulation. If a different ip g results, it must 
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be removed by fusing it with the unpaired parafermion 
mode being moved. The different values used for the 
phases when moving unpaired modes, such as a, b and c 
here, determine how the measurement results are inter¬ 
preted in terms of ip anyons, and so determine the net ip g 
fused with the modes being moved. As such, differences 
in the conventions used for an exchange will change the 
resulting operation only by a factor of T ff , for some value 
of g that depends on a, b, and c. 

For the standard convention used throughout this pa¬ 
per, with a = b = c = 2 , the phase assigned to a ip g 
occupation by the exchange is u> 9 / 2 u} 9 ^ 9+1 ' > . These are 
indeed all square roots of w 9 , as predicted in the main 
text. However, they are not of the elegant form uj 9 / 2 that 
would be more conducive for the proofs of Appendix [D] 

For an exchange operation that does have the required 
form, consider a = 1 and b = c = 2. The phase assigned 
to ip g is w -9 "/ 2 in this case, up to a global phase of w 1 / 2 . 
The required phases uj 9 / 2 would therefore be obtained 
from an anticlockwise exchange. 


The fact that we obtain the phase uj ® 2//2 for a clock¬ 
wise exchange in this case means we would get uo~ 9 for 
a full clockwise nronodromy. This would seem to con¬ 
tradict the arguments of the main text, which predict 
2 

a phase of to 9 . However, note that these phases differ 
only by a factor of ui 2g = w 29 , and so are equivalent up 
to a factor of T 2 . Since such factors are to be expected 
for different choices of a, b and c, this monodromy does 
not contradict our expectations. The ambiguity in these 
phases reflects the fact that the anyon model described 
by the fusion rules in Eq. (141 obeys only projective non- 
Abelian statistics. 

Note that the arguments above do not assume anything 
about the initial state of the exchanged parafermions. 
Only their initial positions and the operations used to 
move them are required. The effect of the braiding is 
expressed in terms of T, the parity operator for their 
shared ip mode, which can be defined for any pair of 
parafermions. It therefore does not matter what state 
the fusion space of the parafermions was initially in, and 
it does not matter whether or not they exist at the end 
of the same defect line. The effect of the braiding is the 
same in all cases. 


Appendix D: Generators of the Clifford group 

For a tensor product of n d-level systems (C d )® n , the 
Pauli group Vd is defined as the group generated by the 
generalized Pauli operators X ( and Z.j , and the Clifford 
group Cd is defined as the normalizer of Vd in the unitary 
group on (C d )®". That is, elements in Cd map tensor 
products of d-level Pauli operators to other such tensor 
products under conjugation. 

We start with some general remarks on the action 
of Cf n on Vf n for d = 4. For d = 4, an opera¬ 
tor X a Z b has eigenvalues {1} if a = b = 0, {1,-1} if 


both a and b are even and at least one of them is non¬ 
zero, {z 1 / 2 ,z 3 / 2 ,z 5 / 2 ,z 7 / 2 } if both a and b are odd, and 
{1, z, — 1, —*} if a + b is odd. Since the number of distinct 
eigenvalues is preserved under conjugation, this implies 
that the Pauli group Vf n decays into distinct orbits when 
Cf n acts on it by conjugation. This is in stark contrast 
to the case where d is an odd prime, which is studied in 
Ref. (771 . where there is only one non-trivial orbit. The 
orbit containing the elements Xi, Z\, ..., X n , Z n consists 
of elements of the form 


u} k+p / 2 Z^X^ ... Z^X b 


(Dl) 


where fc, ai, b\, ..., a n , b n £ Z 4 , at least one of the ex¬ 
ponents ai, bi, ..., a n , b n is odd, and p = "YV-i a V 
determines whether integer or half-integer powers of uj 
appear as phases (we sometimes write uj for i to avoid 
confusion with indices). 

The following proof is an adaption of the proof in Ap¬ 
pendix A of Ref. E 2 , where it is shown that a certain 
set of gates generate the Clifford group Cf n for the case 
where d is an odd prime. The general structure of our 
proof is identical to the one in Ref. while the gener¬ 
ating set and individual lemmas and their proofs are dif¬ 
ferent. After completion of this work, we became aware 
of the more general proof in Ref. mi¬ 
llet us define the single-qudit unitaries 

H =\ii ^ !•?> <■ k \ (° 2 ) 

j,k= 0 


and 


S = | j) 0'| , 


j =0 


T = ^ e ~ i7v/4 , 


1 

2 


j,k 

=0 



Vi 

1 

-Vi 

1 \ 
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Vi 

1 

-Vi 

-Vi 

1 

Vi 

1 

1 

-Vi 

1 

Vi ) 


(D3) 


(D4) 


Lemma 1. The gates , T\ Z\ X, , and \fiH can 
all be generated from S, T, and Z. 

Proof. As S 8 = T 8 = 1, we have S t = S 7 and = 
T 7 . We have V~iH = STS = TST and X = W ZH = 
{V~iHYZ(V~iH). Finally, A' 4 = Z 4 = 1, so X* = X 3 
and Z^ = Z 3 . □ 

If for some Clifford gate U £ Cf n we have 

U{Z^X bl ... Z“"A^)C/ t = aZ^Xi 1 ... Zn n Xn n , 

(D5) 
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with \a\ = 1 , we write 


Let us define a more general controlled operator as 


M(U)(a ll a n ,b n ) T = (ai, b[,... , a' n , b' n ) T . (D6) 

The matrices M(U) £ 2 2nx2ra form a representation of 
Cf n , as M(UV) = M(U)M(V). 

Lemma 2. The gates S, T, and Z generate the entire 
single-qudit Clifford group C® 1 . 

Proof. For some U £ Cf 1 , let M{U) = {%%)■ Preserv¬ 
ing the commutation relations of the single-qudit Pauli 
operators requires that ad — be = 1 (mod 4). One veri¬ 
fies that there are only 48 matrices M in the matrix ring 
Z 2x2 satisfying the requirement detM = 1 (mod 4). We 
have SXS ^ = \JZjXZ and TZT t = y/ojZXf such that 
M(S) = (ii) an( i M(T ) = Once can verify 

by brute force that products of at most 9 factors M(S ) 
and M(T) generate all of the aforementioned 48 matri¬ 
ces. Finally, since XZX^ = wZ and ZXZ t = u>X, we 
can generate arbitrary phases compatible with Eq. (Dl) 
(for n = 1 ). □ 


We define the controlled Pauli-operators 


c st = sf st (c x nc z y. 


It acts by conjugation as 


(D13) 


Z\ r-}c s t Z\ 

X 1 ^c, t ^^XxX^Zl 
Z 2 r^C et 

X 2 ^ Cat Z[x 2 . (D14) 

Up to the phase w st / 2 , this action is identical to the one of 
the conditional Pauli gate Cx s z* studied in Ref. [77. We 
point out again that such a phase is unavoidable for Z 4 , 
as there is, for instance, no unitary U such that Xi 1 — 
X\X 2 Z 2 , since these two operators are not isospectral. 

Let us define the SWAP gate S via S \j) \k) = |fc) | j). 
Evidently, it acts as 


Xi 1 — >s X 2 , Z\ v->s Z 2 , X 2 1 — >s X\ , Z 2 y->s Z\. 

(D15) 


Cx=J2\j) 01 ®* J ' 


l=o 


(D7) 


The gate S thus allows to generate non-local entangling 
gates from nearest-neighbor ones. 


and 

3 3 

Cz = E \i) 01 ® & = E w3k b'> 01 ® |*> (k\ ■ (D8) 

1=0 j,k—0 

Note that Cz has been called A in the main part of this 
work. We write A 1 —B as a shorthand for UAU^ = B. 
We have 

Z\ H>c x Z x 
X\ y^ Cx X x X 2 

y^c x z\z 2 

X-2 y+cx X 2 (D9) 

and 

Z\ y+Cz 
X\ y^c z X\Z 2 
Z-2 y^Cz %2 

X 2 y+Cz Z 1 X 2 , (DIO) 

showing that Cx,Cz £ Cf 2 . 

Lemma 3. The gate Cx can be generated from S , T , 
and Cz- 

Proof. We note that 

HXH t = Z and HZW = AT* . (Dll) 

Thus 

C x = H\C Z H 2 = (ViH 2 )^C z (ViH 2 ) , (D12) 

which together with Lemma [I] completes the proof. □ 


Lemma 4. The gate iS can be generated from S, T, and 

C z - 

Note that the gates S and iS act identically by conju¬ 
gation. 

Proof. Let 

3 3 

Cx ( 1 , 2 ) = E u> 01 ® , C x( 2 , 1 ) = E xi ® li) 01 ■ 

1=0 1=0 

(D16) 

One verifies that 

Cx^C^CxaviViH,) 2 = iS , (D17) 

which together with Lemmas[l]and[3]completes the proof. 

□ 

Let 

P = a P Z^X b f...Z^X b n " 

Q = a Q Z?X*'...Z%'Xfr. (D18) 

All arithmetics involving the exponents aj , Cj, and 
dj that follow are to be understood modulo 4. It follows 
from the commutation relation ZX = ojXZ that PQ = 
w (^.Q)qP where 

n 

{,P i ^ ^ O'idi biCi . 

i= 1 


(D19) 
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Lemma 5. Given P,Q £ Vf n with ( P,Q) = 1, we can 
generate U £ Cf n from S, T, and nearest-neighbor Cz 
such that 


Lemma 7. For any V £ Cf n we can construct U from 
S, T, Z, and nearest-neighbor Cz such that UX\U' = 
VX x V^ andUZiW = VZ 1 V^. 


P^u a P Z<X b 'i ... Z<X b '<’ Pr °°f- Clearly, 

Q^u a Q Z<X< ...Z^X< , (D20) (VZ 1 V f ,VX 1 V^) = (Z 1 ,X 1 ) = 1, (D25) 


with |ap| = |ag| = 1, and there exists j £ {l,...,n} 


such that atdj — 6'c^ = 1. 


Proof. Let P and Q be as in Eq. (D18). Since 


so by Lemma | 6 j we can assume that VXiV^ = X g P' 
and VZiV^ = Z®Q', up to gates that can be constructed 
from S , T, and nearest-neighbor C z . 

Now let 


r. aidi ~ biCi = 1 (D21) 

i=1 

by assumption, there exists j such that 

Tj = ajdj — bjCj £ {+1, —1} • (D22) 

If r : j = 1, we are done. If Tj = — 1, then there is k ^ j 
with rfc = 2 or rfc = — 1. From Lemma [2] we know that 
from Si and we can generate single-qudit unitaries that 
change () in arbitrary ways as long as Ti = aidi — biCi 
is preserved. So up to gates that can be generated from 
Sj and Tj, we can assume that a ; = 0, bj = Cj = 1, and 
dj = 0. If Vk = 2 then, up to gates that can be generated 
from Sk and 7*,, we can assume that a*, = 1, bk = Ck = 0, 
and dk = 2. Finally, if r k = —1 then, up to gates that can 
be generated from Sk and Tk , we can assume that a*, = 1 , 
bk = 1, Cfc = —1, and dk = 2. We note that application 
of a phase-gate Cz(j,k) changes rj to r' = rj + (bkdj — 
bjdk ), and recall that non-local phase gates C Z (j,k) can be 
generated from nearest-neighbor ones and SWAP gates, 
which we can generate according to Lemma [4] In both 
cases (r/c = 2 and r*, = — 1 ), we find that application of 
a phase-gate Cz(j,k ) gives r'- = 1 . □ 

Lemma 6. Given P,Q £ Vf n with (P,Q ) = 1, we can 
generate U £ Cf n from S, T, and nearest-neighbor Cz 
such that 


P i—>■[/ Z P' 


and 


Q^uX®Q ’, (D23) 


with P', Q' £ Vf 


Proof. Let P and Q be as in Eq. (D18). By Lemma[5j we 
can assume that there is j with ajdj — bjCj = 1. Employ¬ 
ing Lemma [4j we can perform a SWAP between qudits 
1 and j. Finally, we perform a single-qudit unitary L 
on qudit 1 which is such that M(L) = ( ~a-)- As 
det M(L) = 1 (mod 4), such a unitary L can be con¬ 
structed from S and T according to Lemma [2] We note 
that 


P’ = a P Z?x!? ... Z^X*" 

Q' = a Q Z?X*‘...Z2'Xfr. (D26) 

We define 

Cst,i = Sf st (c x(hi) nc z{l!i) y , (D27) 

with i £ { 2,..., n}. The gate 

U = (ViHJUQtViHiyUp , (D28) 

with 

n n 

Up = Y[C Kan!i and U Q = Y[C dnCnti , (D29) 

i=2 z=2 


can be constructed from S, T, and nearest-neighbor Cz 
according to Lemmas [T] [3|and [4] 

Using Eqs. ( Dll[ ) and ( |D14 ), we find the sequences of 
mappings 


X\ t-^Up X g P r Z t g P' 

h>c/ q z~ 1 -^=^ a ' d '- biCi) gP' 

X 1+ ^=^ aidi ~ biCi '> (g) P’ , (D30) 

and 

Z\ e-tf/p Zl -^1 ^ U Q X g O' Z®Q', 

(D31) 


up to phases. Using again that XZX' [ = ojZ and 
ZXZ 1 = ujX, and that X can ge generated according 
to Lemma [TJ allows us to generate arbitrary phases com¬ 
patible with Eq. (Dl). Since 


1 = (ZuXi) = (VZ ± V^VX^) = (ZgQ'.XgP') 

n 

= 1 + ' S ^2(a i d i - bid ), (D32) 

i =2 


we finally conclude that 


Z° :i X b > a z Z and Z Cj X dj \-+ L a x X , 

(D24) 


= VX^ 

Z x Z®Q' = VZiV^ , (D33) 


with |ax| = \olz\ = U which completes the proof. 


□ as required. 


□. 
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Theorem 1. Any Clifford gate V £ Cf n can be con¬ 
structed from S, T, Z, and nearest-neighbor Cz- 

Proof. The proof is done by induction over n. The case 
n = 1 is given by Lennna [2] For n > 1, let U be as 
in Lemma pi Since U'V commutes with X\ and Z\, we 
have WV = 1 ® V', where V' £ Cf " -1 acts on qudits 
{2,..., n}. Assuming that the induction hypothesis holds 
for n — 1, V' and hence V can be constructed from S , T, 
and nearest-neighbor Cz- □ 

Appendix E: Defect lines and holes 

Quantum computation in surface codes often uses the 
concept of ‘hole’ defects These are extended ar¬ 

eas in which a single anyon can reside. Their large size 
makes it difficult to measure their anyon occupancy, and 
also to change it without leaving a trace nearby. This al¬ 
lows them to store an additional logical qubit in a topo¬ 
logically protected manner. The code distance is given 
by the size of the hole (for Z errors) and the distance 
to the nearest to its neighbour (for X errors), and so 
can be made arbitrarily large. They are primarily con¬ 
sidered in systems without a background Hamiltonian, 
where they are created and moved using measurements 
[55H55] . However, they can also be created by adiabatic 
means when a Hamiltonian is present [791 - I81] . We now 
discuss this in detail for our system. 

1. Enlarging and shrinking holes 

The stabilizer is generated by the plaquette operators 
M p and E p for all triangular and hexagonal plaquettes. 
Let us consider, however, removing some of these oper¬ 
ators from the stabilizer. In terms of the Hamiltonian, 
this means removing their corresponding terms. 

Specifically, let us remove the plaquette operators E p 
and E q for two triangular plaquettes p and q. This will 
open up a new fourfold degeneracy in the stabilizer space. 
Corresponding Z basis states | g) can be labelled by the 
oj 9 eigenstates of E p . These states are therefore distin¬ 
guished by the type of e g anyon residing in the hole. Due 
to conservation of anyons, the antiparticle e _ 9 must re¬ 
side in q. A further fourfold degeneracy will arise from 
each additional triangular plaquette removed from the 
plaquette. This is because only one removed plaquette, 
such as q, needs to have its occupation determined by the 
conservation of anyons. 

For a logical qudit encoded in the additional stabilizer 
space, an X type operation corresponds to creating a 
particle/antiparticle pair of e anyons. One is placed on 
p and the other on q. The number of sites on which this 
process has support, and so the number of sites on which 
noise must act in order to cause a logical X error, is 
the distance between the two plaquettes. The qudit will 
therefore be topologically protected against such errors 
as long as the plaquettes are well separated. 


The logical Z of the logical qubit corresponds exactly 
to the operator E p , and E q corresponds to Z'. Since 
these are three-body operators, this type of logical error 
requires action on only three qubit pairs. The stored qu¬ 
dit is therefore clearly not topologically protected against 
Z type errors. 

To address this problem, we can deform the lattice 
by making the plaquettes p and q larger. By making 
them arbitrarily large, logical Z errors can be arbitrarily 
suppressed. 

To enlarge p, consider a neighbouring triangular pla¬ 
quette p'. Let us use j to denote the single site shared 
by these. If the stored qudit holds an arbitrary state | g), 
the state of the code will be a uj 9 eigenstate of E p . It will 
also be an uj 9 eigenstate of E p E p f, since the state is a +1 
eigenstate of the stabilizer E p i. 

Consider the adiabatic introduction of the term Xj + 
Xj to the Hamiltonian. This term should be much 
stronger than the adjacent plaquette operators, and so 
will effectively force j into an eigenstate of X and remove 
it from the code. Since this term does not commute with 
E p or E p /, the resulting state will not be an eigenstate of 
these operators. However, the term does commute with 
the product E p E p f since the support of the two plaque¬ 
tte operators on j cancels in this product. It therefore 
remains the same uj 9 eigenstate as it was for the initial 
state. The qudit state | < 7 ) has therefore been effectively 
transferred from the single plaquette p to the combined 
plaquette pp'. The operator E p E p / becomes the logical 
Z, and has support on four qubit pairs rather than three. 
This is further extended as more plaquettes are added us¬ 
ing more Xj + Xj terms. As long as this is done for both 
p and q , the qudit will become topologically protected 
against Z errors as well as X. 

The process used to extend holes can be reversed in 
order to shrink them. Consider a set of triangular pla¬ 
quettes p , p', p" , ... that have been combined into a 
single hole. The basis state | g) of the qubit stored in this 
hole is associated with the uj 9 eigenstate of E p E p 'E p ff ... 
We wish to shrink this hole so that p is no longer a part 
of it. The state will then have a +1 eigenvalue for E p , 
and the qudit state | g) will be associated with the w 9 
eigenstate of E p f E p n ... 

To achieve this, recall that the combination of the pla¬ 
quettes in the hole is is enforced by the strong Xj + Xj 
terms on their shared sites. To remove p from the hole, 
the term E p + E^ should be added to the Hamilto¬ 
nian, and the Xj + X'- term incident on p should be 
removed. Doing this adiabatically will result in a final 
state with the E p term in its ground state, which is its 
+1 eigenspace. Due to conservation of anyon charge, the 
e g anyon that was held in the larger hole must still be 
held in the smaller hole. The qudit state therefore re¬ 
mains |g). 

As well as this being true for each basis state | g), we 
must also be sure that the process preserves coherent 
superpositions. Any process that causes decoherence in 
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this basis will correspond to measurement (by the envi¬ 
ronment) in the Z basis. Any unitary that introduces 
unwanted relative phases can be expressed as a sum of 
powers of Z. As such, these processes must have support 
on all sites on which Z has support, which are all sites 
around the hole. Since the shrinking (and expansion) of 
holes does not have such support, it cannot cause any 
decoherence. 

Corresponding processes can also be applied to hexago¬ 
nal plaquettes. In that case, a logical qudit can be stored 
in the m g occupations of plaquettes. Such qudits can be 
topologically protected by using lines of Zj + Zj terms 
to combine neighbouring horizontal plaquettes. 

Using the processes of extending and shrinking holes, 
it is possible to move them. Gates can then be imple¬ 
mented through braiding. Braiding an e-type hole of tri¬ 
angular plaquettes in state | g) around an to- type hole of 
hexagonal ones in state | h) corresponds to braiding the 
e g anyon held by the former around the m^ of the latter, 
yielding a phase oj gh . This is a qudit generalization of 
the controlled phase gate. 


2. Fusing holes into defect lines 


By considering the alternative stabilizer generators S 
and R discussed in Sec. IV holes can also be created 
which hold ip g anyons. These are formed by similarly 
combining the double plaquettes. Indeed, these are ex¬ 
actly the defect lines considered in the bulk of this paper. 
These have the property that anyons crossing the defect 
line undergo an automorphism that preserves the struc¬ 
ture of the underlying Abelian state: it maps e anyons to 
their dual, the to anyons, and vice versa. This property 
is not shared by the e- and m-type holes. In these cases, 
the lines form a boundary along which one type of anyon 
can condense, but the other cannot cross. It is this differ¬ 
ence that gives the ^-holes additional properties, namely 
the localized parafermion modes at their endpoints, that 
the e- and m-holes do not possess. The topological de¬ 
generacy and protection, however, is a property shared 
by all three. 

An e-type hole and an m-type hole together correspond 
to a two-qudit space. However, let us consider the sub¬ 
space spanned by states \g,g). These are such that the 
e-type hole carries an e g anyon whenever the m-type hole 
carries an m g . A single qudit can be stored in this sub¬ 
space 

Since il> g = e g x m g , two holes as described above hold 
a net ipg. Similar fusion can also be applied to holes, as 
we will now show. Specifically an e-type and an m-type 
hole can be combined into a ^-type hole, and a ^-type 
hole can be split into an e-type and m-type one. 

These processes are in fact a simple generalization 
of the hole extension and shrinking processes described 
above. Suppose we have a defect line along double pla¬ 
quettes P, P\ P", ... This stores a qudit whose basis 
states | g) are oj 9 eigenstates of WpWp/Wp/> ... Let us 


now extend this line. However, rather than adding an¬ 
other double plaquette, we instead add a triangular pla- 
quette p. This can be done by adiabatically introducing 
the strong term X 3 + Xj to the Hamiltonian on a site 
shared by p and the triangular part of P. This term 
does not commute with either E p or Wp. However, it 
does commute with their product. The final state will 
then have the qudit basis states defined by the operator 
EpWpWp'Wp" .... Further such processes can be used 
to extend the e-type part of the defect line. Correspond¬ 
ing processes on the hexagonal plaquettes can be used to 
grow an m-type part of the defect line. An illustration is 
given in Fig. |TT[ 



FIG. 14. Red circles denote defects as in Eq. (131, with 
of parafermions (purple) emerging at the ends of the 
line. Blue circles correspond to terms of the form Xj 


a pair 
defect 


+ x 


while green circles correspond to terms of the form Zj + Zj. 
These grow the e- and m-part of the defect line, respectively. 
Once both the blue and green defect lines have been added, 
the red defect line can be removed. The tp g particle initially 
stored in the red defect line has then been split into its e g and 
m g components. 


When both e-type and m-type parts have been added, 
the original ^-type defect line can be removed. This is 
done simply by removing the parity operator terms along 
its length and allowing the Sp terms to again dominate. 
The end result is that the ij) g originally stored in the 
defect line now resides in the e-type hole as an e g and 
the to- type hole as an m g , corresponding to the state 
\g,g) of their individual qudits. As for the shrinking of 
holes, this process does not have sufficient support to 
distinguish between different basis states. The process 
therefore does not decohere any superpositions of these 
states, nor does it assign any relative phases. 

To recombine the two holes into a single defect line, the 
process is simply reversed. This will be straightforward 
if the two holes are in a state of the form \g, g ), since the 
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state of the defect line will simply become | g). However 
some processes, such as mistakes during error correction, 
could result in holes whose states are not of this form. 
This will introduce frustration that will not allow all of 
the involved triangular and hexagonal plaquettes to re¬ 
turn to their ground state after recombination. 

As an example, consider a state of the form | g, h). This 
corresponds to an e g and an and could arise from an 
initial state | g, g ) if an nih-g were added in error to the 
?n-type hole, or from \h,h) with an e g -h error on the 
e-type hole. 


The anyons e g and an rrif, can combine either to a 
ifjg and rrih-g, or a iph and e g -h- The former will be 
energetically favourable, due to the weaker strength of 
the M p terms. The adiabatic process will therefore re¬ 
sult in ip g being stored on the defect line and an nih- g 
anyon present as an excitation on one of the triangular 
plaquettes that was once part of the holes. Syndrome 
measurement will then detect the m g -h■ However, since 
the value g — h gives information only about the error 
that occurred, and not the value of g or h, this does not 
extract any information about the stored qudit. 



